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Abstract 

Continuous control over azimuthal flow and shear in the edge of the Large Plasma Device (LAPD) 
has been achieved using a biasable limiter which has allowed a careful study of the effect of flow 
shear on pressure-gradient-driven turbulence and transport in LAPD. LAPD rotates spontaneously 
in the ion diamagnetic direction (IDD); positive limiter bias first reduces, then minimizes (produc- 
ing a near-zero shear state), and finally reverses the flow into the electron diamagnetic direction 
(EDD). Degradation of particle confinement is observed in the minimum shearing state and re- 
duction in turbulent particle flux is observed with increasing shearing in both flow directions. 
Near-complete suppression of turbulent particle flux is observed for shearing rates comparable to 
the turbulent autocorrelation rate measured in the minimum shear state. Turbulent flux suppres- 
sion is dominated by amplitude reduction in low-frequency (< 10kHz) density fluctuations. An 
increase in fluctuations for the highest shearing states is observed with the emergence of a coherent 
mode which does not lead to net particle transport. The variations of density fluctuations are fit 
well with power-laws and compare favorably to simple models of shear suppression of transport. 



While flow shear does provide a source of free energy for instability and turbulence, it 
can lead to stabilization of pressure-gradient-driven instabilities and a reduction of turbulent 
transport in magnetized plasmas [HI2]- The transport barrier in the high- confinement mode, 
or H-mode, of tokamak operation |3J is attributed to the spontaneous development of an 
edge flow layer in which strong shearing suppresses transport [U E]. The direct connection 
between the H-mode edge flow layer and improved confinement was first established in 
experiments on the Continuous Current Tokamak (CCT) in which transport barriers were 
generated by directly driving edge flow using torque due to radial currents driven by biased 
electrodes [I]. Biasing has been used to produce improved confinement in a number of 
subsequent experiments including toroidal devices [5H7] and linear magnetized plasmas [8]- 

While ample evidence for transport reduction in the presence of sheared flow exists [TT|[T2] 
and significant effort and progress has been made in developing a theoretical understanding 
of the interaction between sheared flow and turbulence, there are still a number of open 
questions that can be answered by experiment. In particular, the exact mechanism behind 
turbulence modification and transport suppression by shear is still subject to debate: theories 
present a number of mechanisms including radial decorrelation [13], nonlinear reduction of 
turbulent amplitude [14] . and modification of turbulent cross-phase [15] . Evidence for all of 
these mechanisms exists in experimental data [12] , but a comprehensive experimental dataset 
establishing in detail the parameter regimes where each mechanism is important has not been 
acquired. In part, this is due to the fact that most datasets on flow-turbulence interaction 
come from studies of spontaneously generated flow or in cases where precise external control 
over flow and flow shear is not possible. A number of basic plasma experiments have utilized 
biasing techniques to drive flow and flow shear to study flow driven instabilities (e.g. [16[ITT]); 
however, experiments have not been done in which precise external control over flow shear 
has been achieved in higher-density plasmas with drift-wave turbulence to systematically 
study the changes in turbulence characteristics and transport. 

In this letter, we report on the first experiments in which external control of flow is 
used to document the response of turbulence and transport to a continuous variation of 
flow shear, including a zero shear state and a reversal of the flow direction. Shearing rates 
(7 S = dVe/dr, where Vq = E r /B) from zero to up to five times the turbulent autocorrelation 
rate measured at zero flow shear (t~ c 1 ) are achieved. Turbulent particle flux is reduced with 
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increasing shearing rate, regardless of the direction of the flow or sign of the flow shear, 
with significant reduction occuring for 7 S ~ r^ 1 . The observed reduction in particle flux 
is dominated by a decrease in low-frequency (/ < 10kHz) density fluctuation amplitude. 
For low frequency fluctuations, the crossphase between density and azimuthal electric field 
fluctuations remain near zero for all shearing rates. With higher shear (j s > r^ 1 ) we 
observe the emergence of a coherent mode localized spatially in the region of strong flow. 
Fluctuations from this mode appear to increase density fluctuations above 10kHz, but do 
not appear to contribute to particle flux. 

The Large Plasma Device [TS] (LAPD) is a 17m long, ~60cm diameter cylindrical plasma 
produced by a barium-oxide coated nickel cathode. In the experiments reported here, a 
plasma of density ~2 x 10 12 cm -3 and peak temperature of 8eV is produced in a uniform 
solenoidal magnetic field of 1000G. Measurements of electron density, electron tempera- 
ture, and potential (both plasma potential and floating potential) are made using Langmuir 
probes. Measurements of ion saturation current (J sat cx n e y/T^) and floating potential (Vf) 
are taken with a 9-tip Langmuir probe (flush-mount tantalum tips) while temperature and 
plasma potential are determined using a swept Langmuir probe. 7 sa t fluctuations are taken 
as a proxy for density fluctuations for the measurements reported in this work. Density 
profiles are determined by scaling averaged J sa t profiles to line-averaged interferometer mea- 
surements of density. Turbulent particle flux T oc (n e Eg^ is determined through correlating 
density fluctuations from one tip of this probe with azimuthal electric field fluctuations (Eg) 
derived from floating potential fluctuations on two azimuthally separated tips. Azimuthal 
E x B flow is computed using the swept-probe-derived plasma potential. Flows derived 
using this technique compare very well to measurements using Mach probes j9] and flows 
derived from time-delay estimation (TDE) of the velocity of turbulent structures |19j . 

Biasing experiments have been previously conducted on LAPD in which edge profile 
steepening and a reduction in turbulent flux was observed 0, [ID] . In these experiments, 
edge flow was driven by biasing the vacuum chamber wall with respect to the plasma source 
cathode. Transport reduction occurred only for biases above a threshold value. Below the 
threshold, azimuthal flow was localized near the biased wall and no flow or flow shear was 
driven in the region where drift wave turbulence exists. Above the threshold, the flow was 
able to penetrate radially inward; hence, strong flow and flow shear, with shearing rates 
far above the low-flow turbulent autocorrelation rate, was driven in the region of strong 
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density gradient. Recent experiments were successful in achieving more continuous control 
of potential and cross-field flow in the shadow of a small biased obstacle inserted into the 
LAPD core plasma [20]. Both confinement improvement and degradation (formation of 
strong density depletions) were observed with the density profile created by the obstacle in 
this case. 

Motivated by the success of biasing obstacles to control flow, a large annular aluminum 
limiter was installed in LAPD. The limiter provides a parallel boundary condition for the 
edge plasma and is biased relative to the cathode of the plasma source to control plasma 
potential and cross-field flow. The limiter is an iris-like design with four radially movable 
plates located 2.5m from the cathode as shown schematically in Fig. [IJa). The limiters 
create a 52cm diameter aperture; downstream of the limiter, plasma on field lines with 
radial location r > 26cm has the limiter as a conducting end parallel boundary condition 
and plasma on field lines for r < 26cm has the anode/cathode of the source region as a 
parallel boundary condition. An electrically floating conducting end mesh terminates the 
plasma on the far end of the device. A capacitor bank and transistor switch supply a 
voltage pulse to the limiter. The bias pulse lasts 5ms during the flat-top of the ~15ms 
plasma discharge. The limiter is biased from ~10V below to 50V above the anode potential. 
Typically, plasma potential in the core LAPD plasma (plasma on field lines that connect 
to the source region) is very close to the anode voltage and the cathode sits near ground 
(vacuum chamber wall). The anode potential is above the cathode potential by the discharge 
voltage, which was ~40V during these experiments. 

Spontaneous azimuthal rotation of the LAPD plasma is observed when the limiters are 
unbiased (here the limiters are observed to float to a potential ~ 10V below the anode). 
In this state, an edge flow (peaked just outside the limiter edge) is observed in the ion 
diamagnetic drift direction (IDD), as shown in Figure [l](a). Biasing the limiter positively 
with respect to the cathode tends to drive flow in the electron diamagnetic drift direction 
(EDD). As the limiter bias is increased, the flow in the IDD is first reduced, then brought 
to separate near-zero flow and zero flow-shear states, and ultimately reversed with strong 
EDD flow. 

Measurements of profiles of density and particle flux were made for each bias flow state. 
Values are averaged over a range from r = 27cm to r = 31cm, a region where average flow 
and flow shear scale nearly linearly with limiter bias, as shown in Figure [lib). All other 
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FIG. 1: (a) Diagram of the LAPD device showing annular limiter. (b) Velocity profiles using 
plasma potential from swept measurements, (c) Flow at the limiter edge (black, triangles) and 
mean shearing rate, averaged over 27 < r < 31cm (red, circles). 

spatially-averaged quantities shown in this paper are averaged over the same region in space. 



: 

4- 
i 

-Si 
: a - 


* ~ 0.6 

' o 

* c - 0.3 


k! 1 I ■ ' 1 1 ' 1 - 

V i Lim Rel. Anode (V) 
" X »-10V 
StjX »-5V 

i limiter ^fcCT"^" - *-*-, 

edge " 




CO 
CD 
> 

< 




26 31 
Radius (cm) ■ 


-10 





10 20 



Limiter Bias Rel. to Anode (Volts) 

FIG. 2: Density gradient length scale versus limiter bias. Inset shows density profile at three bias 
values. 

Figure [2] shows the variation in the spatially-averaged density gradient length scale, L n = 
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iVlnnl" 1 with increasing limiter bias. As the limiter bias is increased, reducing the IDD 
flow, an increase in the gradient scale length is observed, indicating a degradation of radial 
particle confinement. The gradient scale length peaks when the averaged shearing rate is 
near zero. As the bias is increased further, reversing the flow and again increasing the 
shearing rate, the gradient gradually steepens and the scale length is lowered, indicating 
improved radial particle confinement. 




FIG. 3: (a)Gradient scale length versus shearing rate. (b)Particle flux normalized to no-shear flux 
as a function of normalized shearing rate. Filled symbols represent points with flow in the IDD. 
Inset: Measured turbulent particle flux versus gradient scale length. 



The observed variation of (L n ) with bias is best organized when compared to the shearing 
rate, 7 S , as is shown in Figure |3](a). The shearing rate is normalized to the autocorrelation 
rate of density fluctuations measured in the zero-shear state. An autocorrelation rate of 
t^ 1 rs 28kHz (r ac ~ 36[is) is calculated by taking the half-width at half-maximum of a 
Hilbert transform of the 7 sat autocorrelation function. Confinement improvement (decreased 
(L n )) occurs continuously and gradually with increasing 7,, and reaches saturation for 7^ m 
r^ 1 (a normalized 7 S of 1). The profile steepening appears to be largely independent of the 
direction of the flow (or radial electric field): IDD (filled points) and EDD (open points) 
flow cases follow the same trend when plotted against normalized shearing rate. 

Measured changes in turbulence and turbulent particle flux are consistent with the ob- 
served changes in the density profile. The turbulent particle flux can be written[21j: 
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where n(f) and Eg(f) are the Fourier transforms of the density and azimuthal electric field 
fluctuations; "i{ n ,E e ) is the coherency between density and electric field; and 4>(n,E e ) is the 
cross-phase angle between density and electric field. 



Figure |3^b) shows the spatially-averaged turbulent particle flux as a function of normal- 

ized shearing rate. The turbulent flux decreases continuously with increasing shearing rate; 
however the observed decrease is slightly slower than that observed for L n . The inset in 
Figure |3](b) shows that the variation in turbulent flux is correlated with the changes in L n 
(but scales in a way that is inconsistent with Fick's law using a fixed diffusion coefficient). 
The trend in reduced particle flux is the same for either direction of flow (IDD or EDD). The 
cause for the reduction in turbulent particle flux can be explored by considering individual 
terms in the integrand of Eqn. [TJ 

Density fluctuations were reduced significantly with increasing shearing in these exper- 
iments. Figure |4|a) shows changes in the spatially-averaged density fluctuation spectrum 
with shearing rate. The shearing rate is signed in this figure, and negative shearing rates 
occur for flow in the IDD. Most of the power is located in frequencies < 10kHz and in this 
range, power decreases overall with increasing shearing rate. A decrease of about one order 
of magnitude in fluctuation power is seen between the minimum shear state and the high 
shear regime where L n and particle flux are minimized; this is made clearer in Figure mh). 




At higher shearing rates, 7 S > r^ 1 , a coherent mode emerges. The frequency of the mode 
increases with shearing rate and the fluctuation amplitude is localized to the peak of the 
azimuthal flow. 

Figure [5|a) shows the reduction in total density fluctuation amplitude with shear in two 
frequency bands: all frequencies below 100kHz in black and all frequencies above 10kHz in 
red. With the emergence of the coherent mode, the high frequency fluctuation amplitude 
does show an increasing trend at higher shearing rates but there is a strong overall decrease 
in fluctuation amplitude with shearing. A reduction is also seen in Eg fluctuation ampli- 
tude, as shown in Figure [5](b) ; however this reduction is weaker than observed in density 
fluctuations. The cross-phase between n and Eg does not change significantly with shear- 
ing. As shown in Figure |5j cos[0( n) £ e )] ~ 1 for all shearing rates. For higher frequencies 
(/ > 10kHz), the cross-phase does change with shearing, with cos[0( nj £ e )] trending toward 
zero at higher shear. This crossphase change explains why the coherent mode that emerges 
at higher shearing rate does not contribute to an increase in the particle flux. The coherency 
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FIG. 4: (a) Contour plot of log I sa t /density fluctuation power versus shearing rate and frequency, 
(b) Power spectra for four different values of shearing rate. 

between n and Eg also decreases with shearing rate, as shown in Figure |5j Overall, the de- 
crease in flux is primarily due to a decrease in turbulent amplitude. This observation is 
distinct from previous work with flows driven by vacuum-chamber- wall biasing on LAPD. 
In those experiments, turbulent amplitude decreased little while the turbulent cross-phase 
experienced a significant change, leading to reduced particle flux [10] . In the experiments 
reported here, the magnetic field is higher (1000G versus 400G) and normalized shearing 
rates are lower (near unity). Cross-phase change is expected in cases with very strong shear- 
ing (7 S ^> r^ 1 ) |22j. Future experiments will explore the variation of the turbulent response 
to higher normalized shearing through changing plasma parameters, in particular magnetic 
field. 

Lastly, we add a comparison of our data to a simple theory, the Biglari-Diamond- Terry 
(BDT) model [13], which predicts a power-law scaling with shearing rate of the turbulent 
amplitude of the form: (7 s /r ac 1 ) a . As seen in Figure [5J a best fit of a = 0.530 compares 
favorably to the BDT prediction of a = 2/3 for the reduction in density fluctuation ampli- 
tude. It should be noted, however, that the BDT model is fairly simple and the validity of its 
assumptions is questionable for the experimental conditions reported here. In particular, as 
the shearing rate is increased in LAPD, the density profile is changing (in BDT a fixed drive 



8 



hi i i i i i i i I i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
• • • 0.3kHz-100kHz (yx J" 053 



1.0 



0.5 



0.0 





(b) 8E /6E (y =0) 



i i i i 



-1 

1.0 
0.5 
0.0 



I i i i i i i i i i I i i i 1 1 i i i 



(c) 005(6 







n,E' 

I I 



I I 



(d) n-E Coherency 



' i i i i 



2 3 

(Vs) T ac 



FIG. 5: Components of particle flux versus shearing rate including I sat /density fluctuation 
power (a), electric field fluctuation power (b), crossphase(c) and coherency(d) with black points 
for low or all frequency, red for high only. 

is considered). Future work will focus on direct comparisons to more comprehensive models 
of shear suppression, including comparisons to two-fluid simulations using the BOUT++ 
3D turbulence code [23] . 

This letter presents the first experiments in which the response of pressure-gradient-driven 
turbulence to a continuous variation of shearing rate, including a near-zero flow shear state 
and a reversal in the direction of flow, is studied. Increased shearing improves radial particle 
confinement regardless of the direction of the azimuthal flow or sign of the flow shear. The 
observed reduction of turbulent particle flux with shear is attributed to a reduction in the 
amplitude of density fluctuations. These experiments were performed at a fixed set of plasma 
parameters (fixed magnetic field, neutral pressure, discharge power); future work will explore 
the variation in turbulent response to shear as these parameters are varied. 
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